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Introduction 



For a star, its mass is the most important parameter regarding its evolution 
and interaction with environment. On larger scales, galaxy formation and evolution 
are also shaped by the mass distribution of stars. Moreover the initial star mass 
distribution, the Initial Mass Function (IMF), appears strikingly to be universal 
with a power law behaviour for masses above 1 Mq ( Kroupa||2002 ) and a lognormal 
shape below ( Chabrier et al.pOOO ). Therefore, the IMF is a key issue and a compul- 
sory theory test in star formation, one of the most basic but still most challenging 
problems in astrophysics. 

It is now well established that star formation occurs in the densest and coldest 



parts of molecular clouds, known as dense cores (see Ward-Thompson et al. (1994) 



Kirk et al. (2005), Ward- Thompson et al. (2007)). As cores set the initial conditions 
for star formation, the core mass function (CMF) is often compared to IMF to 
provide constraints on star formation efficiency, timescale and fragmentation. In 
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addition to core density and temperature structures, the CMF may also be a theory 
test by comparing observed and expected CMFs. But above all, the CMF could 
answer the still unsolved question of fragmentation, especially regarding massive 
star formation. The presumed differences between CMFs derived for example from 



competitive accretion model (Bonnell et al. 1997), monolithic collapse (McKee & 



Tan 2003 ) or some in-between models could indeed distinguish these theories. If a 



strong link between CMF and IMF is found, the IMF origin issue may be answered 
by directly solving the CMF one. 

Recent CMF measurements in low-mass star forming regions have shown a great 



likeness in shape between IMF and CMF (Motte et al. (1998) and Johnstone et al. 



(2000) in p Ophiuchi main cloud, Alves et al. (2007) and Rathborne et al. (2009) 



in Pipe Nebula, Testi & Sargent (1998) in Serpens core, [Motte et al^ ( |200lD and 



Johnstone et al. (2000) in Orion B, Tothill et al. (2002) in the Lagoon nebula) 



These observations confirm the role played by cores as direct precursors for stars. 
But it also strengthens the idea of a one-to-one relationship between individual stars 
and individual pre-stellar cores. As a consequence, it enhanced also the idea that 
the shift between the break or peak masses observed in IMF and CMF is due to 
a star formation efficiency less than 1. The star formation efficiency is the ratio 
between the final star mass and its direct progenitor core initial mass; the currently 



favoured value is about 1/3 (Alves et al. 2007). 



The interpretation of the shift between the IMF and CMF as caused by the 



efficiency of star formation is being questioned, however. Swift & Williams (2008) 



have shown through numerical experiments that the comparison between the mass 
functions of dense cores and stars might only distinguish with great difficulties one 
among several evolutionary schemes, given the current observational accuracy. They 
add that the best way of distinguishing several schemes is to study the mass func- 
tions over a wide range of mass scales including the high- or low-mass tails. More- 
over Goodwin et al. (2008) have argued that the CMF derived in Nutter & Ward- 



Thompson (2007) will best match the IMF derived in Kroupa (2002) if cores are 



assumed to form binaries or higher-order multiples. Observational results in Perseus 
dHatchell fc Fuller||2008D have provided some new insights into the possible evolution 
of a prestellar core. The mass distribution of prestellar and protostellar cores in this 
region appeared indeed inconsistent with a simple direct one-to-one mapping to the 
IMF. 

In addition to the nature of the shift between IMF and CMF, the similarity in 
shape between the two distributions could also be interpreted as a consequence of 
the central limit theorem. It has been argued that the IMF is set through the act 
of a large number of independent physical processes and variables, which lead to its 



lognormal shape (Larson (1973), Zinnecker (1984) and Adams & Fatuzzo (1996)). 



Therefore, the CMF would also present a lognormal shape, as it is produced in 
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a similar way. Moreover, the process of deriving a CMF from observational data 
adds independent random processes which differ from those used to obtain the IMF 
(these include random errors on core boundaries, core masses, cores blending which 



do not occur when objects are point sources as in IMF derivation) (see Reid et al. 



2010). The lognormal shape of both IMF and CMF would then be independent of 



a one-to-one relationship between prestellar cores and stars. 

Those uncertainties about the exact relationship between the IMF and CMF 
appear to remain unresolved in a large part because the mass distribution of stars 
that will actually form eventually from an observed CMF is unavailable. In young 
stellar groups, however, where stars have not moved much from their birth sites, 
it is possible to approach numerically what natal cores looked like. The original 
dense core column density map can be estimated based on the current young stellar 
population, allowing the derived CMF to be directly compared to the actual local 
star distribution. 



Starting from YSO distributions divided in groups by Kirk & Myers (2010) in 
four nearby star-forming regions presenting different degrees of crowding, we recre- 
ated a possible original column density map for each group. Each YSO contributes 
to the global column density map of the group by adding one prestellar core profile 
centred at the YSO position. This profile encloses a total mass equal to the YSO 
mass, divided by a fixed star formation efficiency, either 0.3 or 1.0. We used two 
well-known profiles representative of different points of view on star formation, the 



Critically Stable Bonnor Ebert sphere (CSBE) (see Bonnor (1956) and Ebert (1955)) 



and the Thermal Non-Thermal sphere (TNT) (see Myers & Fuller (1992)). Mim- 



icking the observational CMF derivation, we used the clumpfind algorithm clfind2d 



(Williams et al. 1994) to obtain the core mass function and compare it directly to 
the local star mass distribution. From now on, we will address by 'initial' CMF, the 
core mass distribution obtained from the local star mass distribution by a one-to- 
one relationship between cores and stars and by a fixed star formation efficiency; we 
will address by 'derived' CMF, the mass distribution of cores identified by clfind2d 
on the simulated column density maps. The properties obtained from the 'derived' 
CMF will be also called 'derived' (e.g. 'derived' star formation efficiency). 

Our main results are that the numbers of derived cores are undercounted and 
their masses overestimated, due to blending between cores, even in relatively un- 
crowded regions. This happens with no added background and noise and no vari- 
ation in the number of stars per initial core; the inclusion of these effects would 
increase the inconsistency between original and derived cores. We find also that the 
blending effects are so strong in crowded regions that they dominate other effects, 
such as differing core models used to derive the column density maps. 

The data and method used are presented in more details in Section [1} then 
results are displayed in Section [2] and discussed in Section [3] 
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1. Method 



1.1. Initial YSO distributions 



This study is based on the 14 groups derived by Kirk & Myers (2010) from 



catalogs of YSOs in 4 nearby regions. The regions were divided into groups by using 



the Minimum Spanning Tree (MST) algorithm (Barrow et al. 1985), following the 



procedure of Gutermuth et al. (2009 ). The four analyzed nearby star-forming regions 



are Taurus, Chamaeleon I, Lupus 3 and IC 348 and all have catalogs extending to 
late M (or even LO) spectral types. The primary catalog used for each region is 



Luhman et al. (2010) (for Taurus), Comeron (2008) (for Lupus 3), Luhman (2007) 



(for Chamelaeon I) and Lada et al. (2006) plus Muench et al. (2007) (for IC 348). 



For the full catalogs, see Kirk & Myers (2010). 



The distances adopted to each region are the same as in Kirk & Myers (2010) 



140 pc for Taurus, 200 pc for Lupus 3, 160 pc for Cha I and 300 pc for IC 348. 



As stated in Kirk & Myers (2010), these four regions are the only ones within 300 



pc containing stars younger than several million years and whose local extinction 
is low enough to make them observable in the optical/near-IR and have a spectral 
classification completeness of greater than 90 %. 



The mass estimation of YSOs was based on the spectral type procedure outlined 



m 



Luhman et al. (2003) using a combination of models from Palla & Stabler (1999) 



Baraffe et al. (1998) and Chabrier et al. (2000). The completeness limit is about 



0.02 Mq for all regions. 



The MST algorithm used by Kirk & Myers (2010) links all the YSOs of one 



region in a tree, ie a connected graph, which has the property of having the minimum 
total branch length. By using a cut-off length determined by the total distribution 
of branch lengths, the MST algorithm defines these groups in a way that mimicks 



the eye behaviour. See Kirk & Myers (2010) for further details 



More details on the YSO spectral types, positions and mass estimation can be 



found in Kirk & Myers (2010). 



1.2. Going backward : From YSO to cores 



1.2.1. Core models 



A column density map was derived from the locations and masses of the YSOs 
by using different core models and star formation efficiency values. 



The two models used in this study are the Critically Stable Bonnor-Ebert Sphere 
model (CSBE) (see Bonnor (1956) and Ebert (1955)) and the Thermal-Non Thermal 
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model (TNT model) (see [Myers fc Fuller] ( p92| and |Myers| ( |2010[ )). We use these 
particular models as they are frequently used in the community and each offers a 
different point of view on star formation, in particular on the link between the mass 
enclosed and the clump radius. Their exact expressions can be found in Appendix 

m 

The CSBE model was computed in two ways, which we label the CSBET model 
and the CSBEP model. For the CSBET model, the temperature was set to 16 K for 
all groups and the external pressure was allowed to vary. This temperature value is 
the same as the one measured in the Pipe nebula ( Alves et al.|2007 ). For the CSBEP 
model, the external pressure was set to P/k = 3.0x10^ K.cm~^ for all groups and the 
temperature was allowed to vary. This pressure value is thirty times the pressure 
measured in the Pipe nebula, as the Pipe nebula value was too low to describe 
most star forming regions. Typical star forming regions have higher values such 



as Ophiuchus ((Johnstone et al. 2000)) and Orion B (Johnstone et al. 2001). The 
adopted pressure makes the core column density profile peaks at 30x10^^ cm~^, which 



is consistent with high-resolution observations of nearby cores (see [Alves et al.|[200^ 



Kirk et al. 2006). It is noteworthy that these pressure and temperature choices are 



et al. 


2000 


et al. 


2001 



consistent with the ones derived by assuming a Bonner-Ebert profile in (Johnstone 
p Ophiuchi) (T ^ 10-30K and P/k ^ 10^"^ K.cm'^) and ( [Johnstone 
Orion B) (T ^ 20-40K and 3.10^ < P/k < 3.10^ K.cm"' 



For all 



models, we assumed a mean particle mass, m, of 2.3 times the mass of hydrogen. 

Figure [T] displays for each model the column density profiles versus the projected 
radius for an enclosed mass of 0.5, 1, 3, 5 and 10 Mq. Whereas the profiles of the 
CSBET and the TNT models are similar when the mass enclosed is low, they differ at 
high enclosed masses. From now on, we will describe the column density profiles in 
terms of 'core' component (spiky component) and 'clump' component (the extended 
halo surronding the centrally peaked 'core' component). At low masses, the 'core' 
component of the TNT model dominates and is essentially a singular isothermal 
sphere, similar to the CSBE. At high masses, however, the TNT model still shows 
a central 'core' but also a more extended 'clump'. On the other hand, the CSBET 
model retains the same shape but its peak column density scales inversely with 
mass while its extent scales in proportion to mass. Thus as mass increases, the 
CSBET model becomes more extended than the TNT model, while the CSBET 
peak drops below that of the TNT model at the same radius. The CSBEP model 
however doesn't show any important 'core' component. It retains the same shape 
on the entire mass range as the CSBET does but the column density peak value is 
independent of the mass enclosed and is set by the external pressure. 
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1.2.2. Star formation efficiency and Resolution 



The star formation efficiency (SFE) parameter e is defined as the ratio between 
the star final mass and the initial total core mass M where the star comes from. 



M, = eM 



(1) 



We used star formation efficiency values of 1.0 and 0.3 to determine the parent 
core mass for every star. The value of 1.0 was motivated by its use as a reference, 



whereas 0.3 was chosen to match the SFE value derived in the study of |Alves et al. 



(2007) using dust extinction maps. 



Once the parent core mass was obtained, the column density function was com- 
puted on a spatial grid of step 0.001 pc for each YSO and for each different core 
model separately. This step is small enough to allow a good representation of the 
column density proffies (see Figure [T]). The individual column density proffies of 
the progenitor cores within each stellar group were then summed together, to give 
the global column density map necessary to produce the stars seen today. Finally, 
the maps were smoothed to reproduce observational resolution due to submillimeter 



single-dish observations (e.g., the bolometer survey in Perseus of Enoch et al. 2006) 
or star count extinction (e.g., Lombardi et al. 2006 whose dust extinction maps 



were used in Alves et al. (2007)) with a Gaussian kernel of FWHM 0.5 arcmin and 
1.0 arcmin. 

The resulting maps can be compared directly with dust extinction maps since 
the extinction depends primarily on the column density. In contrast, dust emission 
maps are also sensitive to dust temperature; and molecular line maps depend on 
molecular abundance, velocity and line excitation. 



1.3. Identifying observable cores 



We derive the core mass function that would be observed for each group using 
a standard clumpfinding algorithm. We used the two-dimensional version of the 



automated algorithm clumpfind (clfind2d, Williams et al. (1994)), the same one as 



used in the Pipe nebula study of Alves et al. (2007) and Rathborne et al. (2009). 



The algorithm works by first contouring the data at certain levels set by the user. 
Then it searches for peaks of column density which locate the cores. Afterwards the 
algorithm follows these peaks down to lower column density values, until it reaches 
the threshold set by the user. It does not thus assume any core proffie as does 



for instance the Gaussclumps algorithm (Stutzki & Guesten 1990) by assuming a 
Gaussian shape. 



In more detail, at each iteration, the algorithm finds all the contiguous pixels 
which have a value between the current level and the next down one. It then 
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assigns them to a pre-existing core or a new one, depending if they are connected 
or isolated from any previously identified core. In case of blended column density 
features, a 'friends-of-friends' algorithm is used to distribute the pixels between 
several identified clumps. Eventually, at the final level, a core has to be greater than 



a certain number of pixels, if not, it is rejected. For any further details, see Williams 



et al. (1994) 



Our choice of clfind2d parameters was based on Rathborne et al. (2009), who in 



turn relied on a clfind2d performance study by Kainulainen et al. (2009). Rathborne 



et al. (2009) improved on the previous list of cores compiled by Alves et al. (2007) 



using contours of 2x 1.2 mag with a lowest threshold of 1.2 mag. Alves et al. (2007) 



used clfind2d on the Pipe nebula extinction map with a maximum contour level of 
6.0 mag, which led then to the truncation of the extinction contouring at 6.0 mag 
and to the fusion of multiple well-separated extinction peaks into single-extinction 



features. Thus in our study, as in Rathborne et al. (2009), the threshold was set 
to 1.2 xlO^^cm"^ and the difference between two levels at 2 x 1.2 xlO^^cm"^ (we 
used the conversion 1 mag = lO^^cm^^). This parameter choice was motivated by 



the closer comparison to Rathborne et al. (2009) that it offers. 



One of the most important results of the simulations of Kainulainen et al. (2009 ) 



is that the degree of crowding within a molecular cloud can significantly effect both 
the measured core parameters and the derived CMF in a more important way than 
the parameters selected for the core extraction algorithm do. Therefore we separated 
two cases in our study, isolated and blended groups, according to the criterion used 



in the simulations of Kainulainen et al. (2009). According to this criterion, a group 



is isolated if it displays a value of the ratio f of mean separation to mean diameter 
above 1.0. By separation we refer to the distance measured from the peak position 
of a core to the peak position of its nearest neighbor. The separation and the 
diameter values of the cores are defined from clfind2d output (position of the peak 
and surface for each core). For f < 1, Kainulainen et al. ( |2009 ) showed that the 
mass determination of individual cores is very uncertain and that the derived CMF 



may not represent the underlying mass function. Kainulainen et al. (2009) used, for 



the cores, gaussian profiles, random positions and masses, whereas we are using here 
physical profiles, observed positions and masses of stellar groups whose most massive 



stars tend to be in more crowded regions (Kirk & Myers 2010). These differences 



could influence strongly the measured CMF and increase the value of f under which 
the mass determination becomes uncertain; it raises also the question of the quality 
of the determination of the core profile. This ratio f is approximately 2.0 in the 
Pipe Nebula, as it is a relatively quiescent and sparse region. To avoid confusion, 
we denote f as the 'crowding ratio' from now on. 



In Rathborne et al. (2009), the larger-scale background was subtracted prior 
to core identification. Since there is no background added to our maps, and the 
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column density is inferred only from the present-day stars, there is no need to first 
subtract any background. Without background, it is true that the simulated maps 
departed from the observed dust extinction ones. Nonetheless, it allows us to skip the 
background removal step and to study more precisely the effects of the clumpfinding 
algorithm on the actual matter that will eventually accrete on the protostar. Finally, 
it is noteworthy that any additional background will increase the blending of a given 
region. 



2. Results 

2.1. Column density maps 

Figures [2} [3| |4] and [5] display the column density maps of representative groups 
for each region (Taurus 5 or L1529 for Taurus, IC 348 1 or IC 348-Main for IC 348, 
Lupus 3 1 or Lupus 3-Main for Lupus 3, Chal 2 or Chal-South for Chameleon I). 
In Figure [2| the three main steps followed in this study are summarized : starting 
from a YSO group (left panel), we simulate an initial starless column density map 
(middle panel), from which the clumpfind algorithm clfind2d isolated cores (right 
panel) . 

These column density maps show two distinct structural components and are 



then qualitatively similar to dust extinction maps, such as the maps in Kirk et al. 



(2006), where large-scale structure is seen in extinction maps within which small- 
scale features seen in thermal emission are embedded. The simulated column den- 
sity maps similarly show a 'clump' component, which has a spatial extent typically 
around 0.3 pc (but in an extreme case like IC 348 1, it can reach 1 pc) and which 
originates from the external part of several summed column density profiles. On 
top of the clump component appear several 'core' components, which have a spatial 
extent less than 0.1 pc, matching the central part of a single core column density 
profile. This description is particularly accurate when the initial SFE is set to 0.3, 
which is currently the favoured SFE value. 

Regarding clumpfind output, one can easily see by looking at the maps in 
Figures [2] through [5] that cores are undercounted by clumpfind, although the total 
mass within one region is recovered. As expected, this effect is much stronger in 
crowded regions. Figure ^ allows a quick comparison between two extreme cases '. a 
very isolated one (Taurus 5) and a very crowded one (IC 348 1). Whereas Taurus 5 
has a nearly one-to-one relationship between progenitor cores and stars, the crowded 
IC 348 1 shows a much higher number of stars per clumpfind core. From now on we 
will call the number of stars per core the fragmentation ratio F. In IC 348 1, clfind2d 
is less effective in identifying structures: as can be seen in Figure |2| the cores have 
highly irregular shapes and appear unphysical. 
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Figure |3] shows the changes in the column density map of Lupus 3 1 for the 
CSBET model when the values of resolution and initial SFE change. It shows 
qualitatively that fewer cores are found with a decrease in initial SFE and/or poorer 
resolution. A poorer resolution has a stronger effect however than a decrease in 
initial SFE, as the number of identified cores always decreases significantly with the 
poorer resolution. 

Figures |4] and [5] show side by side the column density maps derived for each of 
the three models for the Cha I - South or Cha I 2 group, assuming a initial SFE 
of 0.3 (Fig. |4]) and 1.0 (Fig. |5]), and a resolution of 1.0' for both. The maps show 
a larger difference when the initial SFE is set to 0.3 and when the radii of massive 
cores are bigger. The difference in blending of sources between the panels can be 
explained by the relation between the enclosed core mass M and the core radius R. 
The CSBET model yields core radii that increase proportionally to the total enclosed 
mass, much more rapidly than for the CSBEP model, where R oc a/M and the TNT 
model, where R oc M (see Appendix g in exact form M oc {AR + '-^), but 
the fit R oc M°-^^ is excellent for the present YSOs masses.). For massive stars, 
the CSBET model yields the widest column density profiles, which in turn causes a 
larger degree of blending of cores. 



2.2. Derived core properties 

Several properties can be derived directly from clumpfind output: the total 
number of cores Nc,tot, the mass of each core Mc, the total mass of these cores 
Mc\tot, the radius of each core R and the crowding ratio f describing whether the 
cores appear isolated (f>l) or blended (f<l). The core radius R is derived from 
the core surface returned by clfind2d by assuming a circular projected shape. By 
counting the number of YSOs one core spatially encloses, we derive the average 
number of stars per core, which we denote as the 'fragmentation ratio' F. We fit a 
power-law to the function Mc(R) and derive the slope (3 in the expression McocR^. 
These bulk properties describing each group are given in Table 1. 



2.2.1. Multiple stars per core 

Before comparing the initial CMFs and the derived CMFs, it is interesting to 
calculate the average number of stars found in each core. Table [T] shows that even 
in isolated regions, the ratio F is never strictly equal to one and increases as the 
ratio f decreases, which confirms the legitimacy of the crowding ratio criterion f. 
Figure [6] shows the average value of the ratio F versus the value of the ratio f , where 
each point plotted is averaged over all three models and both resolutions for the 
same group. In addition to the spatial crowding, the number of cores identified also 
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depends on the peak column density. A low resolution or a flatter column density 
profile (as in the CSBEP model) yields fewer cores identified. 

For a given set of parameters and group, the CSBEP model almost always 
leads to the highest fragmentation ratio value of the three models. This may seem 
surprising as the CSBEP model has the smallest radii for massive cores and as 
such, the least blended maps for a initial SEE of 0.3. The high fragmentation can 
be explained by the fact that the peaks of CSBEP profiles are the broadest - the 
CSBEP model does not lead to a 'core' component which is as important as the 
CSBET model for low masses cores or the TNT model for any core (see Figure [T]). 

The CSBET model, however, has the highest fragmentation when the initial 
SEE is 0.3 and thus the core masses are much bigger. The CSBET column density 
peaks are then smaller and broader (see Figure [T]). On the contrary, the 'core' 
component of the TNT model remains narrow in both low and high mass ranges. 
The broadness of the column density peaks has a larger effect on the fragmentation 
than the initial SFE. 

Thus the blending of sources leads itself to undercounting, all the more impor- 
tant when the column density profiles are not spiky enough to lead to a good core 
identification. Such undercounting appears even in groups where the blending is 
very low (e.g., Taurus region). 



2.2.2. Incomplete mass recovery 

Now that it appears that the common observational CMF derivation using 
clfind2d is not reliable for recovering initial cores in blended regions, we want to 
know how the core mass found by clfind2d relates to the input mass of the YSOs 
within the identified core. Similarly we seek the resulting relation on the relation 
between the initial CMF and the derived CMF. The initial CMF corresponds by 
a one-to-one relationship to the local YSO distribution. We know already that, as 
clfind2d uses a threshold, there will always be some material that won't be assigned 
to any core. Hence the total mass of cores will be less than the total mass of gas 
in the column density map. A poorer resolution strengthens this effect as the maps 
are smoother. 

For each derived core, we compared in Figure [7] the mass assigned by clumpfind 
and the mass it should have, given the YSOs it encloses. In isolated regions (Taurus 
and to a lesser extent Chamaeleon I and Lupus 3), for the CSBET and CSBEP 
models, the most massive derived cores are up to 70 % less massive than they 
should be, while the least massive derived cores lack around 30 % of their mass. In 
the same regions, the least massive TNT derived cores get assigned more material 
than they should get by a factor of 30 % too. The most massive TNT derived cores 
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lack up to 50 % of their mass. The differences between the models can be explained 
by the differences in column density profiles : the TNT model always gives a strong 
core component, whereas the CSBEP model shows no strong core component and 
the CSBET model shows only a strong core component at low masses. In crowded 
regions (IC 348 or Chamaeleon I and Lupus 3 when the initial SFE is set to 0.3), 
the mass from the most massive initial cores tends to be spread out sufficiently for 
the massive cores found by clfind2d to miss 50-75 % of their input mass, with some 
of that material being added to the surrounding low and intermediate mass cores. 

The comparison between the initial CMFs and the derived CMFs is thus compli- 
cated by these effects. If low-mass derived cores are less (respectively more) massive 
than they should, the derived CMFs are broadened (respectively narrowed) in this 
mass range. When the blending is important, the masses assigned to derived cores 
are not related to the YSOs they enclose, particularly for the CSBET model which 
yields the most blended maps. 

2.2.3. Core radii 

In the two last sections, we have shown that the derived CMF depends on the 
input model for relatively isolated regions, as the fragmentation ratio F and the way 
a core mass is related to the masses of the YSOs it encloses change with the input 
model (see Figure [T] column 1 and 2, two top panels corresponding to a resolution of 
0.5'). It is then reasonable to think that one can identify the model used to create a 
map, at least in isolated regions. The best way to achieve this is to study the relation 
between the core radius R and the core mass Mc- Each input model is defined by 
a particular power law relationship (CSBET : Mc oc R, CSBEP : Mc oc and 
TNT : Mcoc R^-^ for a SFE of 1 and given the YSO distributions and M^oc R^-^ for 
a SFE of 0.3 and given the YSO distributions, see Appendix [A]). 

Figure |8] shows derived core radius R versus derived core mass Mc, each column 
corresponding to a region and each line to a set of initial SFE and resolution. 

The power law Mc oc R^ was fit by using a least squares method. It appears 
that the derived power law depends on the degree of crowding, the initial SFE and 
the resolution rather than the initial model power law. Strong blending yields a 
higher /3 exponent, since low mass derived cores get assigned nearby gas which are 
actually from the external part of more extended and massive nearby derived cores. 
More strikingly, the value /3 is in most of the cases around 3, which matches a 
constant density profile. 

It appears then that the blending of cores in crowded regions dominates the 
relation between core radius and mass; the input model relationship has little effect 
and one cannot discriminate among initial core models using derived core properties. 
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We show in this section 2^ that the blending, which occurs with no noise and no 
additional background material, leads to a significant undercounting of cores (from 
a factor of 1.4 in isolated groups up to a factor of 20 in blended groups, see Figure 
whose assigned masses do not relate to the YSO they enclose in case of strong 
blending (see Figure [?]). Moreover, core extraction and mass recovery depends on 
the shape of core profile (see Figure [?]) , but the relationship between the mass M 
and radius R of a derived core is independent of the input core profile and obeys a 
power-law of Mcx:R^ (see Figure [s]). 



2.3. Derived core mass distributions 



We now regard the global scale of a group by comparing the mass function of the 
derived cores, the derived CMF, and the mass function of the initial cores, ie the mass 
function of the original YSOs of the group modulo some star formation efficiency,the 
initial CMF. Usually observationally derived CMFs are compared directly to the 
IMF as the local star distribution is not available. Here, however, we do not have 
to assume anything about the future star distribution. 



Figures |9] through 12 show the original YSO and clfind2d derived cores binned 
mass functions for each region and each set of parameters. As discussed earlier, 
the number of cores identified is always smaller than the number of stars. The 
discrepancy becomes more severe in crowded regions (f<l), especially for a poor 
resolution and a low initial SFE. In the less blended configurations, the derived core 
mass distributions and the initial core mass distributions do not show a noticeable 
shift between their peaks. When the configuration becomes more blended, a shift 
between the peaks of the distributions appears, with the magnitude varying with 



the degree of crowding. For instance, in IC 348 for the CSBET model (Figure 12 



first row), the shift between the peaks of the two distributions is about a factor of 
ten, when the initial SFE is set to 1 and the resolution to 0.5', and increases to a 
factor of 30 when the initial SFE is set to 0.3 and the resolution to 1.0'. 

The most important parameters to describe the IMF are the peak mass, where 
the function reaches its maximum (around 0.1 Mq (see Kroupa 2002 Chabrier 



2003)), the break mass, from which the scale free power law behavior seems to 



hold for the high mass range (around 1 Mq (see |Kroupa||2002t |Chabrier||2003D ) and 



the slope of the high mass tail power law (around 2.3 (see Salpeter 1955 Kroupa 



2002 Chabrier 2003[ )). To compare an observed CMF and the IMF, it is usual to 
compare these parameters for both functions, the star formation efficiency being 
defined as the ratio (see equation [T| between the peak masses or the break masses 
of the two functions. The SFE derived from the comparaison between the derived 



13 



CMF and the initial YSO distribution will be called 'derived' SFE. Here we repro- 
duce this comparison for the relationship between the derived CMF and the local 
star mass distribution, matching the initial CMF with the initial SFE shift. 

To make a power law fit of the Differential Core Mass Function (DCMF, dN / dM - 
p{M) oc M~"), we used the maximum likehood estimate (known as MLE, see Clauset 



et al. (2007) for further details, and also Pineda et al. (2009)). This method avoids 



both the problems raised by analyzing the cumulative function and binning the 



data. As pointed out in Mufioz et al. (2007) and Rosolowsky (2005), the upper mass 



limit can induce high curvature at the high mass end in the cumulative function 
shape. Moreover cumulative functions can show curvature leading to a multiple 
power law fit, even in the case where the underlying CMF can be characterized by 
a single power law (see Li et al.|[2007 ). On the other hand, binning the data induces 



most likely information loss (see Rosolowsky 2005). The MLE consists of fitting the 
following function : 

dN a-1 / M 



d- 



(2) 



dM M Break \M Break, 

where Msreak is the minimum mass value at which power-law behavior holds. N^i is 
the number of cores more massive than Msreak and a is the power-law exponent of 
the distribution. The MLE gives an estimate of the exponent and the approximated 
standard error on it : 



a 



i=l 



M, 



-I -1 



M 



Break 



a — 1 



(3) 



(4) 



The lower cut-off for the power-law region, MBreak, was determined by a Kol- 
mogorov Smirnov test. All of these steps were performed using an algorithm based 
on plfit.py, a python implementation of Adam Ginsburg based on the general al- 
gorithm of [Clauset et al. (2007), (see also online |http : / /tuvalu . santaf e . edu/ 
~aaronc/powerlaws/| and for the python implementation agpy.googlecode.com,). 



Figure [13] shows the fit obtained using the MLE method for Cha I 2 for the TNT 
model, initial SFE 0.3 and resolution 1.0'. It appears that the slope obtained for the 
core mass function is shallower than the one obtained for the star mass function. 
In the MLE method, a finite size bias can be present when the number of cores of 



mass above M^j-eafc is under 50 (Clauset et al. 2007). As some groups don't have 



enough cores to allow a good determination of M Break and a, we divided the groups 
into two categories for each initial SFE and core model: blended groups, which have 
a crowding ratio f < 1 for at least one resolution value and isolated groups for the 
remaining ones. The blended groups are listed in Table [2] for each initial SFE and 
core model. In the isolated groups, there are always more than 50 cores identified. 
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In the blended groups, however, there are fewer than 50 cores above Msreak when 
the resolution is 1', for every model and initial SFE. There the fragmentation ratio 
F is high and the total number of cores very low. 



The results of the power-law MLE fit are displayed on Figure [T6| For each 
input SFE (1.0 and 0.3), the derived SFE values are shown in the upper plot and 
the a values are shown in the lower plot. The derived SFE is estimated in two 
ways : SFEpeafc is the ratio of the peak masses of the stellar and the derived core 
distributions and SFE Break is the ratio of the break masses of the stellar and the 
derived core distributions. Mp^ak is estimated as the position of the maximum in the 
mass distributions using a Gaussian kernel because, as this does not require discrete 
data bins, it more faithfully reproduces the detailed structure of the CMF compared 
to a binned mass function ( Silvermaii||1986 ) . MBreak and a are from the MLE power 



law fit. 



The fit results on Figure 16 show how the peak and slope of the derived 
CMF compare to the peak and slope of the distribution of stars. For a one-to- 
one relation between derived cores and stars, with constant SFE, the peak and 
break masses should follow Mpeafc(stars) = (SFEpeafc) xMpeafe (cores), Mb^eafc (stars) = 
(SFEfcreafc) xMbreafc (cores) and the slopes should be identical. Figures [9|-[l2| however, 
show that crowding and poor resolution generally results in a one-to-many relation 
between cores and stars. 

In isolated regions, the initial SFE is recovered by both SFEpeafc and SFEsreak 
as the difference between the initial SFE and the derived SFEpeafc and SFEsreak are 
only a few percent, with the largest difference at a resolution of 1.0'. The mass at 
the peak of the derived CMF exceeds the mass at the peak of the initial CMF by a 
mean factor 1.0. In blended regions, however, SFEpeak and SFEsreak are not similar 
either to the input values or one another. For a resolution of 0.5', the SFEp^ak value 
is around 0.1 and the SFEp^eafc value is around 0.15/0.25 for the three models and 
both initial SFE input. For a resolution of 1.0', all the models for a initial SFE of 
1.0 have a SFEpeafc value around 0.03 and a SFEpreak around 0.06, while all the 
models for a initial SFE of 0.3 have both SFEpeafc and SFEsreak of around 0.1. The 
mass at the peak of the derived CMF exceeds the mass at the peak of the initial 
CMF by a mean factor 12.1. 

Regarding the slopes, all the derived slopes, both for isolated and blended 



groups, are similar to the value of 2.35 from Salpeter (1955) within errors. In 
comparison to the slope derived in the local star distribution, the slopes are slightly 
better recovered in the blended regions. In blended regions, the stellar slope value 
is around 2.4 and the derived core slopes range from 2 to 2.5. In isolated regions, 
however, the stellar slope value is around 3 and the derived core slopes range from 
2.5 to 2.7. 
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To discuss directly the links between the different mass ranges and their prop- 
erties, the resulting cumulative mass functions E(M), the fraction of cores or stars 



with mass greater than M, are shown in Figure 14 for the isolated groups and in 



Figure [15] for the blended groups. Each row corresponds to a different input model 
and each column to a different set of initial SFE and resolution. Each panel shows 
the cumulative functions with the power-law MLE fit for the derived cores, the local 



star distribution shifted by the initial SFE value and the Chabrier IMF (see Chabrier 



2003). The shift between the local star distributions and the Chabrier IMF shows 



that variations from region to region regarding the IMF must be considered when 



comparing IMF and CMF, as stated in Swift & WiUiams (2008). 



In isolated regions, the fragmentation ratio is very close to one, which allows a 
good recovery of the initial SFE by both methods (SFEpeafc and SFE^reafc)- Since 
the fragmentation ratio is not exactly one, some derived cores are actually several 
initial ones blended together. As the more massive initial cores have a more extended 
profile and are almost always positioned in or near the most clustered part of the 



group (see Kirk & Myers 2010), the blending of initial cores happens more often for 
massive cores. Therefore, the positions of Mpeak and MBreak for the derived core 
distribution are very similar to the ones for the star distribution since only the most 
massive initial cores are concerned by the blending in those isolated regions. The 
slope, however, is very sensitive to the high mass tail and the derived CMF slopes 



are shallower than the initial stellar distribution ones. In Figure [Mj, the comparison 
of the star and derived core cumulative mass functions shows that, for a resolution 
of 0.5' (columns 1 and 3), the two distributions only depart in the very last part of 
the high mass tail, as the fragmentation ratio F is very close to one. For a resolution 
of 1' (columns 2 and 4), the fragmentation ratio F is less close to one and the shift 
between the two distributions happens at a lower mass. This does not change the 
derived SFE, but the slopes of the derived cores and stellar distributions differ more 
than in the better resolution case. 

In blended groups, the blending of initial cores is much more important, which 
is a difficulty when trying to derive the SFE and the slope a. The number of 



derived cores can be very small (see Figure 15, column 2, where only a few cores are 
identified). Even when the number of derived cores is sufficient to avoid the MLE 
fitting finite size bias, the derived SFEs are very different from the initial value 
(for instance, the CSBET model with an initial SFE of 0.3 has a SFEpeafc of 0.10 
and a SFE^reafc of 0.25 for a resolution of 0.5', and 0.09 and 0.12 respectively for a 



resolution of 1.0'). As can be seen in Figure 15 in all panels, the star and derived 
core distributions are separated by a shift that is caused by a fragmentation ratio 
F >1. The shift is wider when the resolution is poor - compare columns 1 and 3 
with a resolution of 0.5' and columns 2 and 4 with a resolution of 1'. This effect 
is directly translated in the derived SFE values which can drop up to 70% when 
the resolution changes from 0.5' to 1.0'. As in isolated groups, the derived slopes 
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are shallower than the initial ones: as for the isolated groups, several initial cores 
are blended together producing only one more massive derived core. This happens 
preferentially for the most massive initial cores as they are larger and tend to be in 
a more cluster part of the group. 

The blending, caused by a spatial crowding and a low resolution, seems then 
to be the crucial parameter when deriving the SFE value. It is well recovered in 



isolated groups (within a few percents of the input value, see Figure 14 and 16), 
but not in blended groups, where the SFE is underestimated and independent of 
the input value (around 0.15 for a resolution of 0.5', around 0.05 for a resolution of 
1', the peak mass of the derived CMF exceeds the peak mass of the initial CMF by 
a mean factor of 12.3 see Figure 15 and 16). Regarding the recovery of the slope 



a, it is shallower than the initial one even in isolated regions but can always be 
fitted by power law whose slopes are similar to Salpeter value (2.35) within error 



for both isolated and blended groups (see Figure 14, 15 and 16). The comparison of 



derived CMFs does not distinguish the input models from one another. In addition 
to causing the differences between derived core masses and radii compared to the 
input models, the blending affects also the properties of groups on a global scale, as 
in no case do we obtain a quantitative recovery of the initial CMF by the derived 
CMF. 



3. Discussion 

Starting from YSO masses and positions in 4 nearby star-forming regions, we 
carried out simulations to estimate the initial starless core column density maps 
with a one-to-one relation between YSOs and cores and with different values of 
SFE. The maps are also smoothed with different resolution values. After running 
the clumpfinding algorithm clfind2d on these maps, we derived the CMF for both 
isolated groups and blended groups. If our procedure were 'perfectly recursive', we 
would recover each and every stellar mass that we started with, modulo the SFE. 
In no case, however, do we obtain such recovery because of the blending. The 
blending arise from spatial crowding of cores but also from the smoothing effect of 
the resolution. The blending appears for all three initial core column density models 
considered. 

The difficulty to recover the initial stellar distribution can not be a consequence 
of the choice of parameters regarding the clumpfind algorithm in a significant way. 
The threshold of 1.2xl0^^cm~^ is low enough to describe well all regions, as it can 
be observed for representative groups for each region on Figures 2, 3, 4 and 5 as 
simulated column density maps and clumpfind outputs show similar boundaries. 
Contouring levels of 2x the threshold describe well the column density landscape 



of even the more crowded regions. Figure 17 shows the contouring levels for Chal 
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3 and IC 348 1 in the case of the TNT model, SFE of 0.3 and resolution of 1.0' 
: the difficulty of determining the peaks of different cores is clearly a consequence 
of the blending rather than of the contouring levels value. Using finer contouring 
levels could lead to a significantly better core identification only if they were of small 
fractions of an Av, which is much smaller than observational uncertainties. 

In both isolated and blended groups, it is always possible to fit the high-mass 
tail of the recovered mass function with a power law whose slope is within 0.2-0.3 
of the Salpeter value of 2.35, the errors being in average of ±0.2. These good fits 
do not discriminate among models. The MLE fitting method requires a minimum 
of 50 cores above the break mass to avoid finite size bias, which was not fulfilled 
in blended regions when the resolution is 1.0' and there is strong blending. Despite 
this, the properties derived in those cases are consistent with the properties derived 
in blended regions when the number of cores is sufficient. The fact that we find a 
Salpeter power-law even in the most blended groups where most of the initial cores 
are merged together or when the original YSO mass function shows a steeper high 



mass tail, could be the consequence of the central limit theorem, as is argued in Reid 



et aL (2010). 



In isolated groups (crowding ratio f >1), our procedure is qualitatively sucess- 
ful, but undercounts the cores by a typical factor of 0.7. Their masses are globally 
recovered, as both derived SFEpeak and SFE Break are within a few percents of the 
initial SFE. The fraction of individual core mass recovered, however, can vary sig- 
nificantly. In case of poor resolution (1'), massive derived cores lack up to 70 % of 
their mass. The recovered mass depends on the core density profile and resolution 
as clfind2d is sensitive to spiky and not extended structure. The reason for qualita- 
tive global recovery is that blending effects tend to be local, only a few cores being 
merged together at a time. 

In blended groups (crowding ratio f <1), this procedure is qualitatively and 
quantitatively unsuccessful, because progenitor cores blend into much larger clumps, 
which cannot be parsed into their initial cores. Reasons for this are first that massive 
cores possess a very extended profile and are almost always positioned in a location 



of higher than average stellar surface density (see Kirk & Myers 2010), which leads 
to the blending of core profiles on a more global scale than in isolated groups. 
Second, clfind2d identifies first the peak of a core and then assigns it an extended 
structure. It leads then to the undercounting when individual peaks are merged 
into the summed column density map. As a poor resolution smooths the peaks, it 
worsens also the undercounting. 

These results suggest that, if stars are represented by progenitor cores with a 
fixed mass ratio, the distribution of apparent core masses derived from a column 
density map is an unreliable estimator of the stellar masses, especially in young 
clusters. The blending effects, as the horizontal shift between IMF and CMF, may 
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be easily mistaken with physical effects (low SFE, fragmentation process, different 
evolutionary timescales). 



The blending issue was already addressed by Hatchell & Fuller (2008). In 



Perseus, the mass distribution of prestellar cores shows an excess of prestellar cores 



at low masses. A possible explanation advanced by Hatchell & Fuller (2008) is the 
selection effect due to the blending. This effect is increased by the fact that all 
clumpfinding algorithms look for peaks and can not then identify a core whose peak 
is merged with profile of other cores. The extent of this effect, in comparison to 
the effect of particular evolutionary schemes, as different evolutionary timescales 
for different cores, remains unknown in Perseus. The difficulty in retrieving the 



evolutionary scheme by comparing IMF and CMF was also addressed by Swift fc 



Williams (2008). This difficulty is increased by the impossibility of comparing an 



observational CMF directly to the IMF of stars that will actually form. The present 
results come from the comparison between the local mass star distribution and the 
derived CMF, i.e. the 'observed' CMF, and it is not possible to distinguish different 
evolutionary scheme, apart from the SFE in the isolated groups. 

The smoothed column density maps are similar to observational dust extinction 
maps (see Kirk et al. 2006). Moreover the maps appear to meet the conclusions 



of Smith et al. (2009), particularly in blended case. This study found by using 



smoothed particle hydrodynamic simulations of massive star-forming clumps that 
most of the mass accreted by the massive stars was originally distributed throughout 
the clump at low densities. The radius values in the maps derived here and the results 
in Smith et al. (2009) are also similar, as very massive cores reach radii around 0.5 



pc and low-mass ones have radii less than half a tenth of parsec. 

Observations, however, can also detect cores which will disperse before they 
form any stars. Thus the observed CMF is a mixture of lower-mass cores which 



make no star, and higher-mass cores which make more than one star (see Hatchell 



& Fuller 2008). In our column density maps, we did not add any background or 



additional low-mass cores that will not form stars. Our maps are then less close to 
observations but show the best possible observable relationship between cores and 
stars. Any extra background or lower-mass cores which make no star will increase 
the local blending and strengthens our main results about the crucial role played by 
the blending in the derivation of any CMF. 

Another limitation of our results is that the stars could have moved from their 
birth sites and then the column density maps are not representative of what the 
actual initial column density maps would have been. The motions of stars in the 
first few Myrs, however, are subvirial (see Andre et al.||2007 Offner et al.||2009 ), i.e., 
relatively slow motions and the configuration of stars should not have changed much 
since their birth. Even so, net motions of stars are more outward than inward, so 
our approach probably underestimates actual blending. 
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The maps are simulated by adding the column density profile of every progenitor 
core and thus by assuming an identical formation moment for every star. A way 
of checking the relevancy of the simulated starless column density maps would be 
to do a SPH simulation from the simulated maps back to a YSO distribution and 
compare it to present YSO distributions. This SPH simulation should take into 
account different star formation moments. It would however require assumptions 
about how the original cores are distributed along the line of sight, as we only know 
the YSO positions in the plane of sky. 

By being able to compare directly the derived CMFs to the local stellar dis- 
tribution, our approach allows good estimates of the undercounting of cores due to 
the blending, as we start with a one-to-one relation between stars and simulated 
progenitor cores. The differences between real star-forming regions and the young 
clusters used here underestimate the effects of blending in actual clusters. 

4. Conclusions 

Starting from YSO masses and positions in 4 nearby star-forming regions, we 
carried out simulations to estimate the initial starless core column density maps. 
After running the clumpfinding algorithm clfind2d on these maps, we find that 
an derived CMF can have the right 'shape' to match IMF but can nonetheless 
undercount the number of star-forming cores and overestimate their masses, due 
to blending of initial cores in crowded regions. Such blending can occur even with 
no additional background material, no noise, a 100% star formation efficiency. Our 
results are : 

• In no case do we obtain a quantitative recovery of the initial CMF by the 
derived CMF. 

• The comparison of derived CMFs does not distinguish the input models from 
one another. 

• The derived high-mass tails can always be fitted by power law whose slopes 
are similar to Salpeter value (2.35) within errors for both isolated and blended 
groups. 

• Even in isolated groups, derived cores are undercounted by a factor of 1.4. In 
blended groups, this factor can be as high as 20. 

• The initial SFE is recovered within a few percents in isolated groups, whereas 
it is not recovered in blended groups, where the initial SFE is underestimated 
and independent of the input SFE (around 0.15 for a resolution of 0.5', around 
0.05 for a resolution of 1'). 
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• The peak mass of the derived CMF exceeds the peak mass of the initial CMF 
by a mean factor of 1.0 in isolated regions and by a mean factor of 12.3 in 
blended regions. 

• Mass recovery depends on shape of core profile. 

• The relationship between the mass M and radius R of a derived core is inde- 
pendent of the input cores profile and obeys a power-law of MocR^. 

These results suggest that, if stars are represented by progenitor cores with a 
fixed mass ratio, the distribution of apparent core masses derived from a column 
density map is an unreliable estimator of the stellar masses, especially in young 
clusters. 

Deriving an accurate CMF in blended regions appears to be a very difficult task. 
Clumpfinding programs which rely on peak identification have difficulties identifying 
the extended part of the core profile in blended maps. Even in more isolated groups, 
the outer part of the cores is not well recovered and we find the relationship between 
derived core masses and radii is independent of the input model relationship. Further 
difficulties are expected to arise in real observations where an extended background 
is removed, as the very extended component of the core could be mistaken for the 
background or even be missed if it falls below the detection limit of the observations. 

We have examined the mapping of known YSOs distributions to their modeled 
observable progenitor cores. The conclusions of this work should be very similar 
what observations would give if the local star distribution was available at the same 
time as the core mass function. We expect the results to underestimate the effects 
of the blending, as real observation must deal with background removal and likely 
a slightly more compact configuration of cores than their present-day locations. 

We would like to thank Jaime Pineda for very interesting discussions about 
the clfind2d clumpfinding algorithm, Michael Rcid for sending us his paper before 
publication and discussions about observational bias in deriving CMF and Patrick 
Hennebelle for his comments. We are also thankful to the Smithsonian Astrophysical 
Observatory for partial support of the visit of Manon Michel. 
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Fig. 2. — An illustration of isolated and clustered initial condensations and the 
effect of blending on core identification. The top row corresponds to Taurus 5 and 
the bottom one to IC 348 1, for the CSBET model using a SFE of 0.3 and a blurring 
of 0.5'. Left : The initial YSO distribution is shown with circles of size increasing 
linearly with YSO mass and colored according to the colorbar shown on left (in 
M0). The initial YSO distribution is also displayed in the middle and right figures 
but with markers of smaller size. Middle : The simulated column density map for 
model CSBET with SFE=0.3 and resolution of 0.5'. The blue colorbar indicates 
column density in units of lO^^cm"^. Right : The resulting cores identified using 
clfind2d. Each core's area is shown in a different color. 
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Fig. 3. — An illustration of how the parameters (SFE and resolution) effects the 
Lupus 3 1 group for the CSBET model. Along each row, each map corresponds 
from left to right to (SFE = 1.0, resolution = 0.5'), (SEE = 1.0, resolution = 1.0'), 
(SFE = 0.3, resolution = 0.5') and (SFE = 0.3, resolution = 1.0'). On each map, 
the initial YSO distribution is shown with marker of size increasing with YSO mass 
and colored according to the colorbar shown on left (in M0). Top : The simulated 
column density map. The blue colorbar indicates the column density in units of 
lO^^cm"^. Bottom: The resulting cores identified using clfind2d. Each core's area 
is shown in a different color. 
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Fig. 4. — An illustration of the effect of using different core models in the Chal 2 
group for a SFE of 0.3 and resolution of 1.0'. Along each row, each map corresponds 
from left to right to the CSBET, CSBEP and TNT models. On each map, the initial 
YSO distribution is shown with marker of size increasing with YSO mass and colored 
according to the colorbar shown on left (in Mq). Top : The simulated column 
density map. The blue colorbar indicates column density in units of lO^^cm"^. 
Bottom: The resulting cores identified using clfind2d. Each core's area is shown in 
a different color. 
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Fig. 5. — An illustration of the effect of using different core model in Clia I 2 group 
for a SFE of 1.0 and resolution of 1.0'. Along each row, each map corresponds re- 
spectively from left to right to CSBET, CSBEP and TNT model. On each map, the 
initial YSO distribution is shown with marker of size increasing with YSO mass and 
colored according to the colorbar shown on left (in M0). Top : The simulated col- 
umn density map. The blue colorbar indicates column density in units of lO^^cm"^. 
Bottom: The resulting cores identified using clfind2d. Each core's area is shown in 
a different color. 
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Fig. 6. — The fragmentation ratio F versus the crowding ratio f. Each color represents a 
region (Red : Taurus, Blue : Chamaeleon, Green : Lupus 3 and Black : IC 348) and each 
group is represented by a circle for the SFE value of 1 and by a cross for the SFE value 
of 0.3. Each point is positioned at the average value over the three models for the two 
resolutions while the bars at each point represent the minimum and maximum values for 
the crowding ratio (horizontally) and the fragmentation ratio (vertically) over the three 
models and two resolutions. Strong crowding (f<<l) results in a large number of input 
cores merging in the map, i.e., a high fragmentation ratio F. 
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Fig. 7. — Each column sho'ws for one region the log of the ratio of the mass derived by 
clumpfind for a core (Mc) and the mass the same core should have to produce the YSOs 
it encloses (My 50) versus log(M(7). Chal and Lupus 3 -were regrouped as their properties 
of crowding are much alike. Along each column, each map corresponds respectively to a 
different choice of parameters (SFE, resolution) and displays the results for each model. 
CSBET and CSBEP shows similar trends which are very different from the TNT results. 
These trends seem to disappear as the crowding increases and the SFE decreases. 
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Fig. 8. — Each column shows for one region the log of the core radius against the log of 
the core mass Mc- Chal and Lupus 3 were regrouped as their properties of crowding are 
much alike. Along each column, each map corresponds respectively to a different choice 
of parameters (SFE, resolution) and displays the results for each model. The solid lines 
represent the expected radii given the initial YSO mass distribution of the group, the 
model and SFE inputs. The black dotted lines corresponds to the fit for each model of 
the power law M oc R^. 



- 32 - 




Fig. 9. — The distributions of masses in the Taurus region. Observed cores are 
shown with a black hne (and shading) while the original stellar distribution is shown 
in green. Where the SFE is 0.3, the red line shows the stellar distribution shifted by 
this efficiency. Each row corresponds respectively from top to bottom to the CSBET 
model, CSBEP model and TNT model. Along each row, each figure corresponds 
respectively from left to right to (SFE = 1.0, resolution = 0.5'), (SFE = 1.0, resolution 
= l'),(SFE = 0.3, resolution = 0.5') and (SFE = 0.3, resolution = 1') 
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Fig. 10. — The distributions of masses in the Chamaeleon I region. Observed cores 
are shown with a black hne (and shading) while the original stellar distribution is 
shown in green. Where the SFE is 0.3, the red line shows the stellar distribution 
shifted by this efficiency. Each row corresponds respectively from top to bottom to 
the CSBET model, CSBEP model and TNT model. Along each row, each figure 
corresponds respectively from left to right to (SFE = 1.0, resolution = 0.5'), (SFE 
= 1.0, resolution = l'),(SFE = 0.3, resolution = 0.5') and (SFE = 0.3, resolution = 

1') 
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Fig. 11. — The distributions of masses in the Lupus 3 region. Observed cores 
are shown •with a black hne (and shading) •while the original stellar distribution is 
sho^wn in green. Where the SFE is 0.3, the red line sho^ws the stellar distribution 
shifted by this efficiency. Each ro^w corresponds respectively from top to bottom to 
the CSBET model, CSBEP model and TNT model. Along each ro^w, each figure 
corresponds respectively from left to right to (SFE = 1.0, resolution = 0.5'), (SFE 
= 1.0, resolution = l'),(SFE = 0.3, resolution = 0.5') and (SFE = 0.3, resolution = 

1') 
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Fig. 12. — The distributions of masses in the IC 348 region. Observed cores are 
shown with a black hne (and shading) while the original stellar distribution is shown 
in green. Where the SFE is 0.3, the red line shows the stellar distribution shifted by 
this efficiency. Each row corresponds respectively from top to bottom to the CSBET 
model, CSBEP model and TNT model. Along each row, each figure corresponds 
respectively from left to right to (SFE = 1.0, resolution = 0.5'), (SFE = 1.0, resolution 
= l'),(SFE = 0.3, resolution = 0.5') and (SFE = 0.3, resolution = 1') 
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Fig. 13. — MLE fitting of tlie stellar (left) and core (right) mass function in Cha I 
for the TNT model, with a SFE 0.3 and a resolution 1'. The values fit for the stellar 
mass function are Msreak = 0.49 Mq and a = 2.56. The values fit for the observed 
core mass function are M Break = 1-46 Mq and a = 2.21. 
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Fig. 16. — Star and Core mass distribiitions properties for an input SFE of 1.0 
(Upper graph) and of 0.3 (Lower graph). For each graph are represented: Upper 
paneh The derived values of the SFE by the method of the peak mass (red circle 
marker) and of the break mass (blue triangle marker); the input value of the SFE 
is indicated by a green solid line. Lower panel: The derived value of the slope a 
derived by MLE fitting in the high mass range (bars), the errors on its derivation (red 
errorbars), the number of cores in the high mass range (Number indicated above each 
bar); the Salpeter value of 2.35 is indicated by a solid green line. For both panels, 
the values are ordered according to the degree of crowding (Isolated/ Blended, see 
Table 2), the input model (CSBET/ CSBEP/ TNT) and if the distribution is the 
initial stellar distribution (*) (for the a values), the distribution for a resolution of 
0.5 arcmin (0.5') or the one for a resolution of 1.0 arcmin (1.0'). 
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Fig. 17.— Contouring levels of 2x(1.2xl02icm-2) for Glial 3 (Chal-North) and IC 
348 1 (IC 348-Main) for TNT model, SFE of 0.3 and resolution of 1.0'. The (xy) 
plane matches the spatial position (pc x pc unit) and the z axis matches the column 
density value (in lO^^cm"^ unit). At each level, the contours of the identified clumps 
are drawn. The center of gravity of each current identified clump is marked as a 
black dot. When two previously identified clumps become contiguous at a lower 
level, they are represented as a global new clump and their centers of gravity are 
linked to the one of the new identified clump by a red line. The blending rather 
than the contouring levels leads to a lack of well defined peaks, necessary to the core 
identification. 
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A. Column density in terms of mass and projected radius 



The initial condensation which produces a protostar is represented here by 
analytic expressions for the column density N as a function of total mass M and 



projected radius b, for the truncated isothermal sphere of Bonnor (1956) and Ebert 



(1955), hereafter BE, and for the 'TNT' density profile described by Myers (2010). 

The density profile of the self-gravitating isothermal sphere is approximated to 
a high degree of accuracy by the expression of Natarajan & Lynden-Bell (1997). 

C 



n = Uq 



D 



where uq is the central maximum density. 



(Al) 



(A2) 



is the dimensionless radius, r is the spherical radius and a is the thermal scale length 



{ATrGmnnoy/'^ 



(A3) 



for velocity dispersion a. Here, G is the gravitational constant and m is the mean 
mass per particle, with C=50, D=48,c^=10 and (P=12. 



Integration of the density in equation AT yields expressions for the mass M 
within radius r, and for the column density N at projected radius b. The mass is 



M = 

where the dimensionless mass is 

/X = C I ^ — Arctan 



(47rmno)i/2G3/2 



D [ ^ — Arctan 



(A4) 



(A5) 



If the sphere is truncated at a fixed radius R, the corresponding dimensionless radius 
is denoted ^max = R/a- If the sphere is critically stable, then C,max = 6.46 and 
equation A5 yields /i = 15.85. This case is denoted a CSBE sphere. 

The column density through a truncated sphere of total mass M within radius 
R is then given by 

N 



4 

a jjLV 



27rmG^M 



where 



C 
7 



Arctan 



(<2 _ 02x1/2 
VSmaa; I 



7 



D 



-Arctan 



where 



/5 



h 
a 



(A6) 

(AT) 
(AS) 
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7^ = + 



(A9) 



and 



6^ = d^ + (3^ (MO) 
Note that in equation |A7[ < /3 < S,max- It is useful to write the scale length in 



terms of the total mass by eliminating from equations X3 and X4, giving 

GM 



(All) 



ie : a oc M if 0" (ie T) is set to a fixed value (CSBET model). 



or 



G \"VAA°-' 



ie : a oc M^-^ if P is set to a fixed value (CSBEP model). 



(A12) 



Another density profile is obtained from the assumption that infall is equally 
likely to stop at any moment and from the requirement that the resulting distribution 
of protostar masses follows the initial mass function (IMF). The profile ressembles a 
superposition of 'core' and 'clump' density profiles. It is also similar to the thermal 



and nonthermal (TNT) model of Myers & Fuller ( 1992 ) and to the two-component 



turbulent core model of McKee & Tan (2003) 



n 



(A13) 



where A=34 pc^cm"^ and B=2700 pc^/^cm"^ (see Myers (2010)). 



The mass M within radius R is obtained by integrating equation A13 



M = inm AR + 



7 J 



(A14) 



Similarly, the column density through a TNT sphere of maximum radius R is 

2A(1 + S2)V2 / D \l/3 



N 



R 



-Arctanis) + 2B (^^^-?^^ rfC(l + C')"'/' (A15) 



where s is a dimensionless variable, < s < +oo, and where the projected 
radius is 

b = R/{l + s^y^^ (A16) 
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Table 1:: Star and Core Properties 



vjrIOUp 






X^dl dill. 


'-^S,tot/ '■^C,tot 


(M /M \ 
(_iVl0/ iVl0j 


Himtl Haeriv 


r 


f5 

i 


Taurusl 


B209 


CSBET 


(1,0.5') 


20/15 


10.64/10.38 


1.0/3.1 


1.33 


2.87 


- 


- 


- 


(1,1') 


-/14 


-/9.98 


1.0/3.8 


1.43 


2.25 


- 


- 


- 


(0.3,0.5') 


-/15 


-/34.80 


1.0/2.4 


1.33 


1.57 


- 


- 


- 


(0.3, r) 


-/13 


-/34.33 


1.0/2.7 


1.54 


1.42 


Taurusl 


B209 


CSBEP 


(1,0.5') 


20/15 


10.64/10.27 


2.0/3.8 


1.33 


2.40 


- 


- 


- 


(1,1') 


-/II 


-/9.76 


2.0/4.6 


1.82 


2.36 


- 


- 


- 


(0.3,0.5') 


-/14 


-/34.90 


2.0/2.9 


1.43 


1.76 


- 


- 


- 


(0.3,1') 


-/12 


-/34.36 


2.0/3.7 


1.67 


1.56 


Taurusl 


B209 


TNT 


(1,0.5') 


20/17 


10.64/10.77 


1.4/3.3 


1.18 


2.80 


- 


- 


- 


(1,1') 


-/15 


-/10.34 


1.4/3.3 


1.33 


2.05 


- 


- 


- 


(0.3,0.5') 


-/16 


-/34.13 


1.6/2.9 


1.25 


1.84 


- 


- 


- 


(0.3, r) 


-/14 


-/33.59 


1.6/3.4 


1.43 


1.50 


Taurus2 


L1495E 


CSBET 


(1,0.5') 


30/27 


15.53/15.10 


1.0/3.1 


1.11 


2.36 


- 


- 


- 


(1,1') 


-/21 


-/14.21 


1.0/3.5 


1.43 


2.00 


- 


- 


- 


(0.3,0.5') 


-/24 


-/49.66 


1.0/2.3 


1.25 


1.44 


- 


- 


- 


(0.3,1') 


-/22 


-/49.12 


1.0/2.7 


1.36 


1.23 


Taurus2 


L1495E 


CSBEP 


(1,0.5') 


30/27 


15.53/14.96 


2.0/3.4 


1.11 


2.06 


- 


- 


- 


(1,1') 


-/17 


-/13.87 


2.0/3.7 


1.76 


2.10 


- 


- 


- 


(0.3,0.5') 


-/25 


-/50.92 


2.0/2.7 


1.20 


1.58 


- 


- 


- 


(0.3, r) 


-/23 


-/50.11 


2.0/3.6 


1.30 


1.33 


Taurus2 


L1495E 


TNT 


(1,0.5') 


30/28 


15.53/15.82 


1.4/3.1 


1.07 


2.28 


- 


- 


- 


(1,1') 


-/22 


-/14.99 


1.4/3.3 


1.36 


2.08 


- 


- 


- 


(0.3,0.5') 


-/28 


-/49.79 


1.6/2.9 


1.07 


1.63 


- 


- 


- 


(0.3,1') 


-/24 


-/49.00 


1.6/3.5 


1.25 


1.39 


Taurus3 


B213 


CSBET 


(1,0.5') 


19/16 


8.08/7.84 


1.0/3.9 


1.19 


3.26 


- 


- 


- 


(1,1') 


-/13 


-/7.32 


1.0/3.9 


1.46 


2.76 


- 


- 


- 


(0.3,0.5') 


-/17 


-/26.41 


1.0/2.4 


1.12 


1.83 


- 


- 


- 


(0.3, r) 


-/14 


-/25.85 


1.0/2.9 


1.36 


.69 


TaurusS 


B213 


CSBEP 


(1,0.5') 


19/15 


8.08/7.72 


2.0/3.8 


1.27 


2.92 








(1,1') 


-/12 


-/7.10 


2.0/3.9 


1.58 


2.69 








(0.3,0.5') 


-/16 


-/26.39 


2.0/2.7 


1.19 


1.91 








(o.3,r) 


-/14 


-/25.81 


2.0/3.3 


1.36 


1.69 


Taurus3 


B213 


TNT 


(1,0.5') 


19/17 


8.08/8.43 


1.4/3.2 


1.12 


3.23 








(1,1') 


-/14 


-/7.96 


1.4/3.6 


1.36 


2.58 








(0.3,0.5') 


-/17 


-/25.76 


1.5/2.8 


1.12 


2.04 








(0.3,1') 


-/14 


-/25.19 


1.5/3.4 


1.36 


1.84 


Taurus4 


L1551 


CSBET 


(1,0.5') 


24/18 


22.69/22.29 


1.0/2.8 


1.33 


2.47 
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Table 1 - Continued 



Group # 


Name 


Model 


Par am. ^ 


Ns,tot/Nc,tot^ 


Ms,tot/Mc,tot^ 
(Mo /Mo) 


Pinit/ Pderiv 


F4 


f5 










-/14 


-/21.66 


1.0/3.2 


1.71 


2. in 








(0 3 n 5'") 


-/15 


-/73.28 


1.0/1.9 


i.6n 


1.45 








(0.3,1') 


-/14 


_/72.62 


1.0/2.4 


1.71 


1.45 


Tti 1 1 1* 1 1 

-LCLUl U.9^± 


L1551 


CSBEP 


(1,0.5') 


24/18 


22 fiq/22 2n 


2.0/3.0 


1.33 


2.23 










-/13 


-/21.55 


2.0/3.8 


1.85 


2.19 








(0.3.0.5') 


-in 


-/74.79 


2.0/2.5 


1.41 


1.45 








(0.3. r) 


-/14 


-/74.06 


2.0/3.2 


1.71 


1.74 


1 1 1 1 G ^ 
XcLLLI Uo^i 


L1551 


TNT 


(1,0.5') 


24/18 


22 fiQ/22 3n 


1.4/3.1 


1.33 


2.59 










-/I 5 
/ -'-"3 


-/21 79 


1 4/3 6 


1 60 


2.32 








(Q 3 5'~) 


-/18 


-/73.32 


2.0/2.7 


1.33 


1.53 








(0.3,1') 


-/16 


-/72.62 


2.n/3.2 


i.5n 


1.56 


-LCLUl. i-X O 


L1529 


CSBET 


(1,0.5') 


14/9 


8 25 /8 04 


l.n/3.9 


1.56 


4.27 








\^^^ ) 


-/7 


-/7 fiQ 


1 n/4 6 


2 nn 


3 68 








(V\ 3 5'") 


-/lO 


-/26.93 


1.0/2.3 


1.40 


2.25 








(0.3,1') 


-/lO 


-/26.44 


1.0/3.0 


1.40 


2.06 


Ts* 1 1 r* 1 1 Q 


L1529 


CSBEP 


(1,0.5') 


14/10 


8 25/7 99 


2.0/3.9 


1.40 


4.02 








(I l'~) 


-/7 


_/7.59 


2.0/4.3 


2.00 


3.34 








ro 3 5'"! 


-/in 


-/27 n3 


2 n/2 7 


1 40 


2.41 








(0.3,1') 




-/26.50 


2.n/3.5 


1.56 


1.88 


-l-CLUl. o t_* 


L1529 


TNT 


(1,0.5') 


14/10 


8 25 /8 50 


1.4/3.3 


i.4n 


4.73 








(1 1'~) 


-/8 


-/8.14 


1.4/4.5 


1.75 


3.38 








fO 3 S'") 


-/in 


_/26.44 


1.6/2.8 


1.40 


2.74 








(0.3,1') 


-/in 


-/26.03 


1.6/3.4 


140 


2.34 


T'ei 1 1 ■pi 1 cfi 


L1536 


CSBET 


(1,0.5') 


2,1/27 


17 67/17 25 


1.0/3.6 


1.15 


2.86 








) 


-/16 


-/16.31 


1.0 /3.6 


1.94 


2.48 








fO 3 5'") 


-/26 


-/57.15 


1.0/2.1 


1.19 


1.42 








(0.3,1') 


_/23 


-/56.39 


1.0/2.8 


1.35 


1.40 


TVl 1 1 T* 1 1 Qfl 

-LCLUl U9U 


L1536 


CSBEP 


(1,0.5') 


31/25 


17 67/17 in 

-LI t\J I 1 -LI .-H_/ 


2.n/3.4 


1.24 


2.64 








(1 l'~) 




-/16.16 


2.n/3.8 


1.94 


2.33 








CO 3 5'") 


-/24 


-/57.98 


2.0/2.7 


1.29 


1.69 








/Tl '\ 1 '"1 

IJ . o . 1 J 


/ zu 


/ / .Ul 


z.u/ o.o 


1 .00 


1 .oy 


Taurus6 


L1536 


TNT 


(1,0.5') 


31/25 


17.67/17.77 


1.4/3.2 


1.24 


3.04 








(1,1') 


-/20 


-/17.03 


1.4/3.7 


1.55 


2.30 








(0.3,0.5') 


-/26 


-/57.11 


1.7/2.9 


1.19 


1.81 








(0.3,1') 


-/23 


-/56.29 


1.7/3.3 


1.35 


1.66 


TaurusT 


L1527 


CSBET 


(1,0.5') 


24/21 


16.13/15.69 


1.0/3.0 


1.14 


2.98 








(1,1') 


-/18 


-/14.98 


1.0/3.3 


1.33 


2.12 








(0.3,0.5') 


-/2n 


-/51.49 


l.n/1.9 


i.2n 


1.39 
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Table 1 - Continued 



Group # 


Name 


Model 


Par am. ^ 


Ns,tot/Nc,tot^ 


Ms,tot/Mc,tot^ 
(Mo /Mo) 


Pinit/ Pderiv 


F4 


f5 








(0.3,1') 


-/19 


-/50.73 


1.0/2.4 


1.26 


1.29 


TVinrii^T 


L1527 


CSBEP 


(1,0.5') 


24/21 


16 13/15 55 


2.0/3.0 


1.14 


2.61 








\^^^ ) 


-/17 


-/1 4 74 


2 0/3 6 

Z(aU / OaU 


1.41 


1.92 










-/21 


_/52.77 


2.0/2.3 


1.14 


1.56 








(0.3,1') 




-/51.83 


2.0/3.0 


1.26 


1.31 


Tanrii^T 

-LCLUX U O 1 


L1527 


TNT 


(1,0.5') 


24/22 


16 13/15 97 


1.4/2.9 


1.09 


2.92 








y^i^ ) 


-/18 


-/15.24 


1.4/3.4 


1.33 


2.14 








(0.3.0.5') 


1 


-/51.96 


1.7/2.5 


1.14 


1.78 








(o.3,r) 


_/20 


-/51.14 


1.7/3.0 


1.20 


1.56 


-LclUl U90 


L1517 


CSBET 


(1,0.5') 


16/15 


12 27/11 94 


1.0/2.4 


1.07 


2.31 










-/13 


-/11.48 


1.0/3.1 


1.23 


1.90 








Co 3 5'") 


-/13 


-/38.58 


1.0/1.7 


1.23 


0.97 








(0.3,1') 


-1X1 


-/38.44 


1.0/2.2 


1.33 


1.01 


1 (1 Li 1 Li o ij 


L1517 


CSBEP 


(1,0.5') 


16/15 


12 27/1 1 90 


2.0/2.5 


1.07 


2.22 










-/13 


-/1 1.43 


2.0/3.5 


1.23 


1.90 








fO 3 5'") 


-/14 


-/40.32 


2.0/2.2 


1.14 


1.41 








(0.3,1') 


-/14 


-/39.84 


2.0/3.0 


1.14 


1.29 




L1517 


TNT 


(1,0.5') 


16/15 


12 27/12 08 


1.4/2.7 


1.07 


2.53 










-/13 


-/II. 59 


1.4/3.4 


1.23 


1.98 








fO 3 5'") 

\\J t\J <^\J tKJ I 


-/15 


-/40.20 


2.0/2.5 


1.07 


1.49 








(0.3,1') 


-/14 


-/39.74 


2.0/3.2 


1.14 


1.37 


Chall 


Chal 


CSBET 


(1,0.5') 


12/10 


3 66 /3 47 


1.0/3.8 


1.20 


2.70 




-SW 




) 


-/8 


_/2.98 


1.0/4.0 


1.50 


1.71 








(0.3.0.5') 


-/II 


-/II. 82 


1.0/2.5 


1.09 


1.43 








(0.3. r) 




-/1 1.44 


1.0/2.8 


1.20 


1.25 


Chall 


Chal 


CSBEP 


(1,0.5') 


12/10 


3 66 /3 39 


2.0/3.9 


1.20 


2.31 




-SW 






-/6 


- 12.11 


2.0/4.2 


2.00 


1.81 








cn 3 n 5'") 


-/in 


-/II. 78 


2.0/2.7 


1.20 


1.44 








(0.3,1') 


-/lO 


-/11.40 


2.0/3.4 


1.20 


1.22 


Chall 


Chal 


TNT 


(1,0.5') 


12/11 


3 66/4 11 


1.4/3.0 


1.09 


2.32 




GUST" 






/Q 


/ O.DO 


1 4 fi 
1.^/ o.O 


1 .oo 


1.(0 








(0.3,0.5') 


-/II 


-/11.93 


1.5/2.4 


1.09 


1.54 








(0.3,1') 


-/lO 


-/11.57 


1.5/3.2 


1.20 


1.30 


ChaI2 


Chal 


CSBET 


(1,0.5') 


96/63 


40.45/39.18 


1.0/3.2 


1.52 


1.30 




-South 




(1,1') 


-/34 


-/37.08 


1.0/3.3 


2.82 


1.22 








(0.3,0.5') 


-/64 


-/131.15 


1.0/2.8 


1.50 


0.66 








(0.3,1') 


-/41 


-/130.86 


1.0/3.2 


2.34 


0.68 


ChaI2 


Chal 


CSBEP 


(1,0.5') 


96/59 


40.45/38.91 


2.0/3.6 


1.63 


1.22 
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Table 1 - Continued 



Group # 


Name 


Model 


Par am. ^ 


Ns,tot/Nc,tot^ 


Ms,tot/Mc,tot^ 

(Mq/Mq) 


Pinit/ Pderiv 


F4 


f5 










-/34 


-/37.08 


2.0/3.8 


2.82 


1.18 








(0 3 n 5'") 


-/56 


-/133.04 


2.0/3.2 


1.71 


0.81 








(0.3,1') 


-/37 


-/131.87 


2.0/3.8 


2.59 


0.83 


ChaI2 


Chal 


TNT 


(1,0.5') 


96/67 


4f) 45 742 38 


1.4/3.2 


1.43 


1.27 










-/44 


-/40.81 


1.4/3.4 


2.18 


1.12 








(0.3.0.5') 


-/69 


-/132.67 


1.6/3.5 


1.39 


0.85 








(0.3. r) 


/47 


/131.41 


1.6/3.7 


2.04 


0.79 


ChaI3 


Chal 


CSBET 


(1,0.5') 


43/30 


21 69721 1 1 


1.0 /3.6 


1.43 


1.39 




-North 






-/18 
/ ^° 


-/20.22 


1 073 5 

± .yj / kj.o 


2 39 


1 16 








fO 3 5'"! 


-In 


-/70.42 


1.0/3.0 


1.59 


0.89 








(0.3,1') 


-121 


-/70.21 


1.0/4.0 


2.05 


0.79 


ChaI3 


Chal 


CSBEP 


(1,0.5') 


43/25 


21 69720 99 


2.0/3.9 


1.72 


1.43 




-North 




(1 V) 


-/13 


-/20.00 


2.0/4.6 


3.31 


1.26 








(V\ 3 5'") 


-/25 


-/71.45 


2.0/3.7 


1.72 


1.03 








(0.3,1') 


-/16 


-/70.83 


2.0/3.9 


2.69 


0.94 


ChaI3 


Chal 


TNT 


(1,0.5') 


43/33 


21 69721 90 


1.4/3.2 


1.30 


1.39 




-North 






-/19 


-/20.96 


1.4/3.3 


2.26 


1.14 








fO 3 5'"! 


-/34 


-/69.99 


1.6/3.3 


1.26 


1.00 








(0.3,1') 


-719 


-/69.19 


1.6/3.7 


2.26 


0.96 


Tjimim^ 


TjiiFiim^ 


CSBET 


(1,0.5') 


36/22 


18 18717 55 


1.0/3.3 


1.64 


1.79 




-]\^ain 




(1 l'~) 


-714 


-/16.37 


1.0/3.5 


2.57 


1.68 








fO 3 S'") 


-/22 


-/59.28 


1.0/2.7 


1.64 


1.03 








(0.3,1') 


-718 


-/58.62 


1.0/3.7 


2.00 


0.99 


T .iim i<£5t 

■Li U^UOtJ 




CSBEP 


(1,0.5') 


36/20 


18 18717 44 


2.0/4.1 


1.80 


1.81 




— IVIairi 




) 


-712 


-/15.91 


2.0/3.6 


3.00 


1.53 








fO 3 5'") 


-/22 


-/59.61 


2.0/3.9 


1.64 


1.23 








(0.3,1') 


-716 


-/58.75 


2.0/4.1 


2.25 


1.17 






TNT 


(1,0.5') 


36/25 


18 18718 44 


1.4/3.1 


1.44 


1.73 




— IVIain 




(1 l'~) 


-715 
7 


-/17.15 


1.4/3.7 


2.40 


1.67 








CO 3 5'") 


-/26 


-/59.07 


1.7/4.1 


1.38 


1.29 








(C\ '\ 1 '"1 
IJ . o . 1 ) 


/I R 

7 -Lo 


/ OO. 10 


1. / / o.y 


9 nn 

Z.UU 


1 1 R 
i.±D 


IC3481 


IC348 


CSBET 


(1,0.5') 


186/41 


111.86/110.59 


1.0/5.2 


4.54 


0.68 




-Main 




(1,1') 


/7 


-/110.02 


1.0/4.9 


26.57 


0.57 








(0.3,0.5') 


/35 


-/370.89 


1.0/5.0 


5.31 


0.48 








(0.3,1') 


-/6 


-/370.51 


1.0/4.2 


31.00 


0.48 


IC3481 


IC348 


CSBEP 


(1,0.5') 


186/32 


111.86/110.63 


2.0/5.2 


5.81 


0.75 




-Main 




(1,1') 


-/6 


-/llO.ll 


2.0/5.7 


31.00 


0.53 








(0.3,0.5') 


-/21 


-/371.93 


2.0/7.0 


8.86 


0.61 
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Table 1 - Continued 



Group # 


Name 


Model 


Param.^ 




Ms,tot/Mc,tot' 
(Mo /Mo) 


Pinit/ Pderiv 


F4 


f5 








(0.3,1') 


-/5 


-/371.36 


2.0/4.7 


37.20 


0.47 




IC348 


TNT 




1 Sfi /44 


111 Sfi/ll 3 Q3 


1 4/5 3 


4.23 


73 




-iVltilll 










1 4/4 S 

± .'4: / 4:.0 


9"^ 9^ 


yJ.'JO 










-/37 




1 7/5 3 


5 03 










(0.3,1') 


-/7 


-/3fi5 30 


1.7/4.2 


26.57 


0.64 




IC348 


CSBET 


(1,0.5') 


11/5 


3 1 3/2 8Q 


1.0/5.2 


2.20 


0.97 




-North 






-/I 




10/49 


1 1 00 










fo 3 n V] 


-/4 


-/10.22 


10/5 


9 75 


80 








ro 3 1 '"i 


-/2 


-/I n on 


10/49 


5 50 


79 




IC348 






ll/o 


o. lo/ Z.oo 


Z.U/ O.Z 


o.D / 


110 
i. iy 




-North 




(1,1') 


-/I 


-/2.60 


2.0/5.7 


11.00 










(0.3,0.5') 


-/3 


-/10.19 


2.0/7.0 


3.67 


0.90 








(0.3,1') 


-/2 


-/9.97 


2.0/4.7 


5.50 


0.76 


IC3482 


IC348 


TNT 


(1,0.5') 


11/4 


3.13/3.34 


1.4/5.3 


2.75 


1.02 




-North 




(1,1') 


-/I 


-/3.08 


1.4/4.8 


11.00 










(0.3,0.5') 


-/5 


-/9.98 


1.5/5.3 


2.20 


0.76 








(0.3,1') 


~/2 


-/9.77 


1.5/4.2 


5.50 


0.79 



Table 2:: Blended groups 



Model SFE 1 SFE 0.3 



CSBET IC 348-Main, IC 348-North L1517, Chal-South, Chal-North, Lupus 3-Main, 

IC 348-Main, IC 348-North 
CSBEP IC 348-Main, IC 348-North Chal-South, Chal-North, IC 348-Main, IC 348-North 
TNT IC 348-Main, IC 348-North Chal-South, Chal-North, IC 348-Main, IC 348-North 



^ (initial star formation efficiency, resolution) 

^ Total number of stars/ total number of derived cores and total mass of stars/ total mass of 
derived cores within each group (stars data from Kirk fc Myers||20T0 1 

^Slope of McxR*^ for the input model/Slope of McxR*^ for the derived cores.. 

^Crowding ratio f = mean separation/mean diameter 

^Fragmentation ratio F = mean number of stars per observed core 



